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We apply a recently proposed theoretical concept and numerical approach to obtain time delays in extreme 
ultraviolet (XUV) photoionization of an electron in a short- or long-range potential. The results of our numerical 
simulations on a space-time grid are compared to those for the well-known Wigner-Smith time delay and 
different methods to obtain the latter time delay are reviewed. We further use our numerical method to 
analyze the effect of a near-infrared streaking field on the time delay obtained in the numerical simulations. 
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1. Introduction 

Recent experimental observations revealed substantial time delays in photoionization of atoms 
from different orbitals 0^3] ■ For some of these observations the so-called attosecond streak 
camera technique 0] has been used to map the desired time information onto the momentum or 
electron spectrum of the photoelectron as a function of the delay between the ionizing extreme 
ultraviolet (XUV) pulse and a streaking near-infrared (NIR) pulse. Theoretical analysis of the 
experimental observations has mainly focused on two aspects recently P, 0, 15r-[l8| . On the one 
hand it is the calculation of the time delay, which is assumed to being related to the so-called 
Wigner- Smith (WS) time delay 5110], which is a well-known concept in scatte ring theory and 



different methods have been proposed for the calculation [H, H, H, 12, 3- On the 

other hand there is the question how and to which extent the observation of the time delay is 
influenced by the NIR field, which is used for the observation in the experiment [H-S HI, 17 , 18], 



In order to contribute to the theoretical analysis we have recently proposed a theoretical 
concept along with a numerical app roach to determine time delays in XUV photoionization 
with and without streaking field [2l|]. In this approach we make use of a fundamental quantum 
mechanical definition for the time a particle spends in a certain region of a potential. The 
definition is analogous to that used in the derivation of the so-called WS time delay [H, H3] . 
Furthermore, we proposed a back-propagation method by which the concept can be applied 
in actual numerical solutions of the time-dependent Schrodinger equation (TDSE) on the grid. 
Application and first results for XUV photoionization in the short-range Yukawa potential and 
the long-range Coulomb potential have been presented recently (2l| . 

In this paper we extend our studies and analyze the XUV photoionization in further short- 
range potentials, such as the Hulthen potential and the Woods-Saxon potential, as well as 
the combination of a short- and a long-range potential. The results allow us to establish the 
applicability of our numerical approach in more general. Furthermore, we use the results of our 
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numerical calculations for a comparison with the results for the WS time delay. This enables 
us to discuss and review different approaches to obtain the WS time delay for short- as well as 
long-range potentials. Finally, we present an analysis of the impact of a NIR streaking field on 
the numerical results for the time delay in different potentials. Our results show, in general, that 
the effect is small as long as the intensity of the NIR field does not exceed 1 x 10 13 W/cm 2 . 



2. Numerical simulations of time delays in XUV photoionization 



In this section we discuss the application of our recently proposed numerical approach [21[ to 
determine time delays in photoionization of an electron by an ultrashort XUV pulse. To this 
end, we first briefly outline the theory as well as the different potentials used in the present 
calculations. We then present results of test calculations showing some characteristic features of 
our time delay calculations for short- and long-range potentials as well as the combination of 
both. 



2.1. Outline of theoretical approach 

We define the time delay in XUV photoionization as the difference between the time t^^R an 
initially bound electron needs to leave a certain region R of a potential due to the interaction 
with the XUV field and the time %(o> r a free particle spends in the same region [2l| (Hartree 
atomic units, e = m = h = 1 are used throughout the paper): 

Ai#,.R = t^fl — t^ ) >R , (1) 

with 

f°° f \\$j( ion )( r t)\ 2 



oo 



and 



t*(o),H= f dt f dr\^°\r,t)\ 2 . (3) 

J~oo JR 

is the initial bound state of the electron in the potential U(r) and is the free state 
of the electron corresponding to the ionizing part of the wavefunction ^( ion ) after transition 
from the initial state to the continuum by absorbing an XUV photon. In Eqs. ([2]) and ([3]) the 
wavefunctions $(* on ) and \E r (°) are renormalized in order to get physically reasonable times. 

The definition given in Eq. J^Q) is based on the well-known concept of a time delay in particle 
scattering of a potential [19l . 1201 ] . As in these early studies in scattering the time delay At^/j 
has a finite limit as the radius of R grows to infinity for short-range potentials U(r). In contrast, 
for long-range potentials such as the Coulomb potential At^^R diverges in this limit 2l(. Fur- 
thermore, according to the definition above the time delays are negative for propagation of the 
particle in an attractive potential U(r), since a free wave packet spends more time in a given 
region R than the corresponding wave packet that has the same asymptotic energy propagating 
in the attractive potential. 

In order to use the above definition in a numerical simulation of an XUV photoionization 
process we use the back-propagation technique [ill ]. To this end, we first solve the corresponding 
TDSE of the system, initially in the state *$>(r,t = 0), under the interaction with the external 
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light field on a space-time grid: 

i-y(r,t)=(^ + V(r) + Vn giA (t))*(r,t), (4) 

where p is the momentum operator and V\i s ht(t) represents the interaction with the ionizing 
light field, while V{r) includes the (short- or long-range) electrostatic potential as well as, if 
used in the calculations, a streaking- field Streak- We proceed with the (forward) propagation 
to large distances on the grid at least until the interaction with the ionizing light field as well 
as the interaction with the streaking field ceases. This allows us to spatially separate the ion- 
izing part of the wavefunction from the remaining bound parts. Alternatively, we could project 
onto analytically or numerically known bound states for the separation. Then we propagate the 
remaining ionizing part of the wavefunction backwards in time without taking account of the 
interaction with the light field V^ight once within the potential V(r): 



dt 

and once as a free particle: 



= (y + W) * ( ' lon) (r,t) (5) 



i^(0)(r,t) = ^*(°)(r,t) (6) 



The corresponding times tig r and t^m R can then be calculated for regions R with radius smaller 
than the location of the wave packet at the start of the back-propagation. 



2.2. Model systems and numerical simulation 

Below we present numerical results for the application of this theoretical method to the photoion- 
ization of an electron initially bound in several different model potentials, namely the long-range 
Coulomb potential in ID: 

V c (x) = - J_ , (7) 
\J x l + a\ 

with Z = 3.0 is the effective nuclear charge, and a\ = 2.0 is the soft-core parameter; short-range 
interactions represented by the ID Hulthen potential Vn(x) and the ID Woods-Saxon potential 
Vws(x), given by (for studies using the short-range Yukawa potential, see [2j|): 

V R , ws (x) = , (8) 

q _ e v , x 2 +a 2 /b 2 

with Vq = 0.5, q = 1.0, 02 = 1.5, and 62 = 6.0 for the Hulthen potential and Vq = 5.0, q = —1.0, 
d2 = 1.85, and 62 = 3.0 for the Woods-Saxon potential, and a combination of Vq(x) and Vws{x)'- 



with xq = 30.0 or 100.0, and b\ = 1.0 for the combined potential. For these potential parameters, 
we obtained almost identical ground state energies for the potentials, namely —1.7118 (Coulomb 
potential Vc(x)), —1.7088 (Hulthen potential Vh(^)), —1.7100 (Woods-Saxon potential Vws(^)) > 
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and —1.7118 (combined Coulomb- Woods-Saxon potential Vc-ws(x))- We note that the combined 
potential has the same ground state energy as the Coulomb potential, since the Woods-Saxon 
factor in Eq. ([9]) represents a smoothed step-function, which does not influence the short-range 
part of the Coulomb potential. 
For the interaction with the XUV light pulse we used length gauge, i.e. 

Vi igU (t) = E XVY (t)x (10) 

where Exxjv represents a linearly polarized pulse with a sin 2 envelope, i.e., 

Exuv(t) = E sin 2 (irt/T)sm(ujt + <p), (11) 

where Eq is the peak amplitude, r is the pulse duration, uj is the central frequency, and 4> 
is the carrier-envelope phase (CEP). The streaking field Kjtreak? if used in the calculations, is 
represented in an analogous form. 

To solve the corresponding TDSE in our simulations, we used the Crank-Nicolson method in 
a grid representation, with, in general, a spatial step of Sx = 0.1 and a time step of St = 0.02. 
In case of the short-range potentials (incl. the combined Coulomb- Woods-Saxon potential), the 
grid extended from —9000 a.u. to 9000 a.u., while a larger grid of —14000 a.u. to 14000 a.u. was 
adopted to forward propagate the ionizing wave packet for the Coulomb potential. The larger 
grid in case of the long-range potential assures that the final momentum of the electron at the 
end of the forward propagation is close to the asymptotic one, which reduces the error in the 



calculation of the time the free particle spends in the potential 211 ]. 

The initial ground states were obtained by imaginary time propagation. For the back prop- 
agation we propagated the ionizing parts at negative and positive x backwards in time inde- 
pendently, either under the influence of the potential, or as a free particle. We determined the 
corresponding times t\f, t R and t^(o> ^ for both parts of the ionizing wavefunction and added the 
two contributions. In the ID calculations we defined the region as R = [±Xi nner , ±cc outer ], where 
dinner and x outcr are the inner and outer boundaries, respectively, and the ±-signs apply to back- 
propagation of the two parts of the ionizing wave packet along the positive/negative x-axis, 
respectively. We absorbed the wavefunction beyond the inner boundary Xi nncr using the exterior 



complex scaling method 22j, |23( ] 



2.3. Some characteristic features 

In Figs. Q] and [2] we present some characteristic results obtained for the time delay At^/j in 
the different potentials without streaking field. As can be seen from the results in Fig. [IJ At^/ t R 
depends on the size of the region it! and is negative in all cases studied. It converges to a 
finite limit as R increases towards infinity for short range potentials onl y. T hese results are in 
agreement with the general features of a time delay in scattering theory [19j, |20[] . In the case of 



the combined Coulomb- Woods-Saxon potential (blue solid and green dashed lines in Fig. Ufa)), 
one can most clearly see that the time delay converges for distances beyond the effective range 
of the potential. Our results for the Coulomb potential (red dash-dotted curve in Fig. QJa)) also 
confirm the expectation that there is no limit for the time delay as the outer boundary increases 
in such a long-range potential. This is due to the well-known logarithmic divergence related 
with ionization from any bound state within the potential. The effect of the long-range tail of 
the Coulomb potential can be clearly seen from the comparison of corresponding results (red 
dash-dotted line) with those for the combined Coulomb- Woods-Saxon potentials in Fig. [2(a). 
We however note that for any finite distance the time delay is a well-defined finite number for 
each of the potentials. This allows us e.g. to study effects of an ultrashort streaking pulse even 



for the long-range Coulomb potential with the present theoretical approach [211 ]. 
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Figure 1, Time delays as a function of the outer integration boundary Pouter for different potentials: (a) Coulomb potential 
(red dash-dotted line) and combined Coulomb- Woods-Saxon potential (blue solid line for xq = 30, green dashed line for 
xq = 100), (b) Hulthcn potential (blue solid line) and Woods-Saxon potential (green dashed line). In all calculations we 
have used an XUV pulse with peak intensity J = 1 X 10 15 W/cm 2 , central frequency uj = 100 eV, pulse duration r = 400 
as, and carrier-envelope phase <j> = for the ionization. (The color version of this figure is included in the online version of 
the journal.) 








Figure 2. Time delays as a function of inner integration boundary Xi nncr . Three short-range potentials are considered here: 
Hulthcn potential (red dash dotted line), Woods-Saxon potential (black dotted line), and the combined Coulomb- Woods- 
Saxon potential (blue solid line for xo = 30 and green dashed line for xq = 100). Laser parameters are the same as in Fig. 
Ill The corresponding WS time delays were shown as black circles. (The color version of this figure is included in the online 
version of the journal.) 



We also note from the results in Fig. Q] that in each case the absolute value of the time delay 
increases most strongly in the region close to the center of the potential, where the potential is 
strongest. Therefore, our theoretical results strongly depend on the choice of the inner boundary 
dinner of the region R, which can be seen from the results in Fig. [2j In the latter calculations we 
fixed the outer boundary of R at x ou ter = 500, which is large enough to obtain converged results 
for all short-range potentials considered here. In our studies presented below we have chosen 
dinner = 0, which corresponds to the expectation value of x for all the bound states considered 
here. As indicated in Fig. [2] and discussed below, for this choice of Xi nncr our results are also in 
excellent agreement with those for the so-called WS time delay, if applicable (black circles in 

Fig. ED. 
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3. Numerical results and discussion 

In this section we first present comparisons of our numerical results for time delays with those 
for the so-called WS time delay for short-range potentials. Then, we use the fact that all time 
delays are finite for any finite region R to study the effect of a NIR laser pulse, as used in recent 
streaking experiments [lj], on the time delays. 



3.1. Comparison with Wigner- Smith time delays 

For particle scattering of a potential the time delay, defined analogous to the expression given 
in Eq. ([1]), can be given in the limit of an infinitely large region R as [lH ]: 

A*w S = %, (12) 

where <p is the phase shift induced by the potential V(r) upon scatterin g o f the particle. This 
expression is also commonly known as the (asymptotic) WS time delay [3, 20] and is well defined 



for any short-range potential. In contrast, application of the concept and formula in case of a 
long-range potential is questionable, since the required limit of the time delay does not exist 
(c.f., e.g., O). 

We therefore restrict ourselves to a comparison of our numerical results for the time delay with 
those for the WS time delay in the case of short-range potentials. We obtained the WS time 
delay for photoionization in two ways. On the one hand we made use of Eq. (I12j) and determined 
the phase derivative in a time-independent scattering approach. To this end, we first considered 
scattering of an electron, incident from x = — oo with a momentum k of the short-range potential 
of interest. We solved the corresponding time-independent Schrbdinger equation numerically 
using the fourth order Runge-Kutta method up to |x| = 500, i.e. well beyond the effective range 
of any short-range potential considered here. Then we projected the numerical solution onto 
the appropriate plane wave solutions for x — > =Loo, and obtained the WS time delay for the 
scattering process At^^ as the derivative of the phase shift in the plane wave propagating 
in positive x-direction with respect to the energy of the incident particle. For a specific XUV 
photoionization process we averaged the results over the energy spectrum of the wave packet, 
as obtained in our time-dependent numerical simulations. Finally, we divided the result by two 
to take account of the fact that photoionization can be considered as a half-scattering process. 

We also employed an alternative concept to determine the WS time delay, which has been used 
in several recent theoretical works P, [l(| • It is based on the assumption that - significantly far 
away from the center of a short-range potential - the trajectory of the electron depends linearly 
on the propagation time. The time delay is then given as the time difference between the time 
obtained by extrapolating the linear part of the electron trajectory back to the center of the 
potential and the time zero. Usually, the electron trajectory is calculated as the expectation value 
of the ionizing wave packet as a function of time. Here, we instead obtained the trajectory by 
plotting the outer boundary x ou tor of the region R as a function of t^.^ (blue solid lines in Fig. 
[3]). In case of a short-range potential the trajectory beyond the effective range of the potential 
was then fitted to a linear line. As an example, we present in Fig. [3]^a) the results for the 
combined Coulomb- Woods-Saxon potential with xq = 30. It can be seen that the extrapolation 
does not depend on the range over which the trajectory is fitted (cf., green dashed and red dash- 
dotted lines) and the WS time delay can be determined from the intersection of the extrapolated 
curves with the base line (see inset of Fig. El^a)). In contrast, this extrapolation method is not 
applicable in the case of Coulomb potential. Due to the logarithmic divergence in the trajectory 
a fit to a linear line is not possible. Linear fits to different parts of the trajectory indeed lead 
to a significant variation in the extrapolated time delay (see green dashed and red dash-dotted 
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Figure 3. Extrapolation method to calculate the WS time delay. The solid blue lines show the outer integration boundary 
2'outcr as a function of time for (a) the combined Coulomb- Woods-Saxon potential with xq = 30 and (b) the Coulomb 
potential. The green dashed and red dash-dotted lines are obtained by linearly fitting the solid blue lines in two regions: 100 
to 200 and 400 to 500 respectively. The two insets show the behaviors near t = for each line. At in panel (a) corresponds 
to the WS time delay for a short-range potential. (The color version of this figure is included in the online version of the 
journal.) 



lines in Fig. EJb)). This is in agreement with the fact that the (asymptotic) WS time delay is 
an ill-defined concept for a long-range potential such as the Coulomb potential. 

Our numerical results are in excellent agreement with those for the WS time delay in case of 
the short-range potentials. For example, in Fig. [2] the WS time delays are shown as black circles 
and are found to agree well with the numerical results for Xi nner = 0. Next, we compare in Figs, 
d] and [5] our numerical results (blue diamonds) for the time delay as a function of XUV frequency 
at a fixed pulse duration of r = 400 as and as a function of the XUV pulse duration at a fixed 
central frequency of oo = 100 eV with the WS time delays obtained either via the calculation of 
the energy derivative of the phase shift (dashed lines) or the extrapolation method (red open 
squares). One may notice that there is a small discrepancy between the time delays from the 
time-dependent simulations and the WS time delay calculated as the phase shift derivative from 
the time-independent approach, which can be most clearly seen in Fig. 03(c). This discrepancy 
is related to our choices of grid parameters. Test calculations have shown that this discrepancy 
can be removed by further reducing the spatial step size. 

The latter results also agree well with qualitative expectations. The absolute value of the time 
delay decreases towards zero as the frequency of the ionizing XUV laser pulse and, hence, the 
final kinetic energy of the emitted electron increase, since the effect of any potential becomes 
more and more negligible (see panels (a)-(d) in Fig. 0]for the different short-range potentials). 
Furthermore the absolute value of the time delay decreases with an increase of the XUV pulse 
duration (see Fig. [5]), which is closely related to the dependence on the XUV frequency. Due to 
the finite pulse duration the time delay obtained for the wave packet can be considered as an 
average over contributions at particular energies within a certain bandwidth (weighted by the 
ionization probability at a given energy). Since the time delay does not change linearly with 
the kinetic energy, the value obtained for a wave packet will be smaller than its contribution at 
the expectation value of the kinetic energy. This difference decreases and, thus, the time delay 
for the wave packet increases as the energy bandwidth of the wave packet decreases, i.e. as the 
pulse duration increases. Besides, the expectation value of the kinetic energy of the ionizing 
wavepacket increases as the XUV pulse duration increases, which also causes the time delay to 
increase. 
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Figure 4. Time delays as a function of XUV central frequency uj: (a) combined Coulomb- Woods-Saxon potential with xo = 
30, (b) combined Coulomb- Woods-Saxon potential with xq = 100, (c) Hulthen potential, and (d) Woods-Saxon potential. 
In each panel, three sets of time delays are shown: TDSE results (blue diamonds), WS delays using the extrapolation 
method (red open squares), and WS delays obtained via the phase shift derivative (dashed lines). Other laser parameters 
are: I = 1 X 10 15 W/cm 2 , r = 400 as, and <j> = 0. 



3.2. Streaking effects 

Finally, we study the effect of an additional external NIR field on the time delay obtained in 
our numerical simulations. This plays an important role in view of recent observations of time 
delays using the attosecond streaking camera technique [H, 0]. In the latter a weak NIR field 
is used to map time onto momentum and retrieve from the momentum (or, energy) spectrum 
of the photoelectron as a function of the delay between ionizing XUV pulse and streaking NIR 
pulse information about time delays between different processes. 

As mentioned above we take account of the streaking field as one of the potential terms in V(r) 
in Eq. @ and therefore consider the streaking field on equal footing with the short- or long-range 
electrostatic potential of interest. In this way we obtain as result of the back-propagation step 
the time delay as 

A (NIR) _ (NIR) , v 

where t^ 1 ^ ^ is the time the ionizing wavepacket spends in region R in the presence of the electro- 
static potential and the NIR streaking field. In our calculations we consider XUV photoionization 
from the ground state of the short-range potentials and checked that the ionization induced by 
the NIR field is negligible up to an streaking field intensity of 1 x 10 13 W/cm 2 in each of the 
cases presented below. 

In order to study the effect of the streaking field, we present in Fig. [6] the difference between 



January 16, 2013 2:26 Journal of Modern Optics time'delay 



9 



-1.1 
-1.2 
7-1.3 
ra -1.4 

0) 
TD 

o> 

■I ~ 1 ' 5 
-1.6 

-1.7* 



-1.5 



TDSE 

phase shift derivative 

□ extrapolate 



(a) 







JB' 


♦ TDSE 




phase shift derivative 


/ 


□ extrapolate 


f 

I 

H 

< 


(b) 



200 400 600 800 1000 1200 1400 1600 



200 400 600 800 1000 1200 1400 1600 



-0.65 
-0.7 
~3 -0.75 
I -0-8 



■I -0.85 
-0.9 



O TDSE 

— phase shift derivative 
□ extrapolation 



(c) 



200 400 600 800 1000 1200 1400 1600 
XUV pulse duration (as) 



-0.9 



-1.2 



-1.3 



♦ TDSE 

phase shift derivative 

□ extrapolation 



(d) 



200 400 600 800 1000 1200 1400 1600 
XUV pulse duration (as) 



Figure 5. Time delays as a function of XUV pulse duration r. All the symbols have the same meanings as in Fig. H Other 
laser parameters were: I = 1 X 10 15 W/cm 2 , uj = 100 eV and <j> = 0. (The color version of this figure is included in the 
online version of the journal.) 
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Figure 6. Time delay differences AT as a function of IR intensity: (a) the combined potential with xq = 30 (blue line 
with diamonds) and xq = 100 (green line with triangles), (b) Hulthen potential (blue line with diamonds) and Woods- 
Saxon potential (green line with triangles). The time delay Aim induced by the IR is shown as red open squares in panel 
(a). The XUV pulse is added at the middle (zero field, maximum vector potential) of the IR field. Laser parameters are: 
7 XUV = 1 x 10 15 W/cm 2 , u xuv = 100 cV, r XU v = 400 as, fcuv = 0, / N IR = 1 x 10 12 W/cm 2 , A NIR = 800 nm, N mK = 3 
cycle, 4>nir = 0. (The color version of this figure is included in the online version of the journal.) 

the time delays in the streaking field to those without streaking field, i.e. 

AT = At^ R) - At*,*. (14) 



We obtained converged time delays for all short-range potentials once the NIR field ceased. 
In each of the present calculations the XUV pulse was applied at the center of the NIR field, 
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which corresponds to the maximum of the vector potential. Our results confirm our previous 



conclusions from studies using the Yukawa and Coulomb potentials [21[ that the time difference 
increases with increase of the intensity of the streaking field, but the effect remains small up to 
streaking field intensities of about 1 x 10 13 W/cm 2 . 

It is however interesting to note that the variation of the time delay with the streaking field 
intensity can be quite different for different short-range potentials. This can be seen, in particular, 
from the results for the combined Coulomb- Wood-Saxon potentials with different effective range, 
shown in Fig. [6|a) . We obtain a considerable smaller increase in the time delay as a function of 
intensity for the potential with a larger effective range (xq = 100, blue line with diamonds) than 
for the potential with shorter effective range (xq = 30, green line with triangles). This can be 
understood as follows. The effect of the NIR streaking field onto the time delay can be thought 
of as adding two contributions to the original delay At^/j: a time delay caused by the IR field 
only (A£nir) and a time delay induced by the coupling of the NIR field and the electrostatic 
potential (AtNiR-Ep), i-e. 

AiJ5f } = At* ii? + At N iR + At N iR_ EP . (15) 
By combining Eq. (|14p and (|15p . we get for the time delay difference: 



AT = AtNiR + Ai N iR_ EP . (16) 

In our numerical simulations we can determine Ai^iR by back-propagating the ionizing part of 
the wavefunction in the NIR field alone, determining the corresponding time and subtracting 
the time for the free particle back-propagation t^m R . Since these test calculations in the back- 
propagation step do not account for the specific electrostatic potential and the ground state 
energies for each of the potentials considered are similar, our results for AiNiR are almost the 
same for each of the potentials (for a given set of XUV and NIR streaking field parameters) and 
an example is shown as red open squares in Fig. [6^a). 

The time delay A£nir induced by the NIR field agrees very well with the full time delay 
difference AT for the combined potential with xq = 30 (blue line with diamonds in Fig. [6^a)), 
which indicates that the effect of a coupling between NIR streaking field and electrostatic po- 
tential is negligible in this case. This makes sense since we applied the XUV ionizing field at the 
maximum of the vector potential (and, hence, a zero of field) and the streaking field does not 
change significantly from zero while the ionizing wavepacket propagates within the electrostatic 
potential, as long as its effective range is small. As the effective range increases (xq = 100) the 
coupling between electrostatic potential and (now varying) streaking field becomes stronger and 
the time delay difference decreases (green line with triangles in Fig. E|a)). 



4. Conclusions 



We have applied a recently introduced numerical approach to obtain time delays in numerical 
simulations to XUV photoionization from different (short- and long-range) potentials. We have 
found, in general, excellent agreement with the results of calculations for the so-called WS 
time delay if the potential is short-ranged. In contrast, our numerical results confirm that for 
a long-range potential, such as the Coulomb potential, the time delay diverges in the limit of 
an infinitely large region and, hence, the WS time delay does not exist. Finally, we have found, 
in general, that the impact of a streaking field on the numerical results for the time delay is 
negligible as long as the intensity of the streaking field does not exceed 1 x 10 13 W/cm 2 , which 
corresponds to the intensity limit at which ionization of the target by the streaking field itself 
sets in. 
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